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NL63 coronavirus (NL63-CoV), a prevalent human respiratory virus, is 
the only group | coronavirus known to use angiotensin-converting 
enzyme 2 (ACE2) as its receptor. Incidentally, ACE2 is also used by 
group II SARS coronavirus (SARS-CoV). We investigated how different 
groups of coronaviruses recognize the same receptor, whereas ho- 
mologous group | coronaviruses recognize different receptors. We 
determined the crystal structure of NL63-CoV spike protein receptor- 
binding domain (RBD) complexed with human ACE2. NL63-CoV RBD 
has a novel f-sandwich core structure consisting of 2 layers of 
B-sheets, presenting 3 discontinuous receptor-binding motifs (RBMs) 
to bind ACE2. NL63-CoV and SARS-CoV have no structural homology 
in RBD cores or RBMs; yet the 2 viruses recognize common ACE2 
regions, largely because of a “virus-binding hotspot” on ACE2. 
Among group | coronaviruses, RBD cores are conserved but RBMs are 
variable, explaining how these viruses recognize different receptors. 
These results provide a structural basis for understanding viral evo- 
lution and virus—receptor interactions. 


receptor protein | SARS coronavirus | 
spike protein receptor-binding domain | virus-binding hotspots 


fundamental yet unresolved puzzle in virology is how 
viruses evolve to recognize their receptor proteins (1). 
Specifically, how do different viruses recognize the same recep- 
tor protein, and how do similar viruses recognize different 
receptor proteins? Do viruses select their receptor proteins by 
chance, or do they target specific virus-binding hotspots on these 
receptor proteins? Structural information of virus—receptor 
interfaces can potentially answer these questions. To date, 
although a few studies have obtained structural information for 
a single virus—receptor interface (2—6), no study has provided 
structural information for the interfaces between different vi- 
ruses and their common receptor protein. Here we provide 
such structural information, by showing that nonhomologous 
receptor-binding proteins of 2 coronaviruses bind to the same 
“virus-binding hotspot” on their common protein receptor. 

A recently identified human coronavirus, NL63 (NL63-CoV), is 
associated with common colds, croup, and other respiratory dis- 
eases (7, 8). Potent neutralizing antibodies against NL63-CoV are 
detected in sera from nearly all humans older than 8 years, 
suggesting that NL63-CoV infection is common in childhood (7, 9). 
NL63-CoV belongs to the coronavirus family, a group of enveloped, 
positive-stranded RNA viruses that infect many mammalian 
and avian species. Coronaviruses are classified into 3 serologic and 
genetic groups: mammalian group I, mammalian group II, and 
avian group III (10). NL63-CoV is the only group I coronavirus 
known to use angiotensin-converting enzyme 2 (ACE2) as its 
receptor (9), whereas the others use aminopeptidase-N (APN) 
(10-12). Curiously, ACE2 is also the receptor for the severe acute 
respiratory syndrome (SARS) coronavirus (SARS-CoV) (13), a 
group IT coronavirus responsible for SARS (14, 15). 

Coronaviruses enter cells through a large spike protein on their 
envelopes (10). The coronavirus spike protein is a membrane- 
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anchored trimer and contains 2 subunits, receptor-binding subunit 
S1 and membrane-fusion subunit S2 (Fig. 14). The S2 subunits 
from group I and group IJ coronaviruses share both sequence and 
structural homology (16); they contain homologous heptad-repeat 
segments that fold into a conserved trimers-of-hairpin structure, 
which is essential for membrane fusion (Fig. 14) (16, 17). Surpris- 
ingly, the S1 subunits from group I and group IJ coronaviruses have 
no obvious sequence homology. Nevertheless, they can both be 
divided approximately into N-terminal region, central region, and 
C-terminal region (Fig. 2A and B). The S1 central regions of both 
NL63-CoV and SARS-CoV are defined receptor-binding domains 
(RBDs) that are sufficient for high-affinity binding to their com- 
mon receptor ACE2 (18-22). 

To date, the crystal structure of SARS-CoV RBD complexed 
with ACE2 is the only atomic structure available for any 
coronavirus S1 (Fig. 2D) (5). SARS-CoV RBD contains 2 
subdomains, a core and a receptor-binding motif (RBM); RBM 
exclusively contacts ACE2. ACE2 contains a claw-like peptidase 
domain, with 2 lobes encircling the active site. SARS-CoV binds 
to the outer surface of the N-terminal lobe of the ACE2 
peptidase domain. Structural information has been lacking for 
group I coronavirus S1, either alone or in complex with its 
receptor. 

How do NL63-CoV and SARS-CoV both use ACE2 as their 
receptor, despite no obvious sequence homology in their S1 
subunits? One hypothesis is that NL63-CoV and SARS-CoV 
share homologous RBMs, and that through RNA recombina- 
tion, SARS-CoV acquired its RBM from NL63-CoV or an 
NL63-CoV-related group I coronavirus, gaining binding affinity 
for ACE2 and infectivity for human cells (20, 23). Extensive 
mutagenesis studies have been done to characterize the inter- 
actions between NL63-CoV and ACE2 but have failed to yield 
consistent results regarding the NL63-CoV binding site on ACE2 
or the domain boundaries of NL63-CoV RBD (19-21). It is also 
intriguing that NL63-CoV and other group I coronaviruses 
recognize different receptors despite obvious sequence homol- 
ogy in their S1 subunits. Here we report the crystal structure of 
NL63-CoV RBD complexed with its human receptor ACE2, 
revealing structural mechanisms whereby NL63-CoV recognizes 
the same receptor as SARS-CoV but a different receptor from 
other group I coronaviruses. 
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Structure of NL63-CoV RBD complexed with human ACE2. (A) Domain structure of the NL63-CoV spike protein. Unique, unique region; NTR, N-terminal 


region; Central, central region; CTR, C-terminal region; HR-N, heptad-repeat N; HR-C, heptad-repeat C; TM, transmembrane anchor; IC, intracellular tail. The 
unique domain only exists in NL63-CoV and is not involved in receptor binding (19, 20). (B) Overall structure of NL63-CoV RBD complexed with human ACE2. The 
RBD core isin cyan, RBMsin red, and ACE2 in green. (C) Averaged electron density map contoured at 1.0 cand covering a portion of the NL63-CoV—ACE2 interface. 
(D) Sequence and secondary structures of NL63-CoV RBD. Beta-strands are drawn as arrows. RBMs are in red; the remainder of the RBD is in cyan. Disordered 
regions are shown as dashed lines. (E) Kinetics and binding affinity of NL63-CoV RBD and human ACE2 by surface plasmon resonance using Biacore. Structural 


illustrations were made using Povscript (31). 


Results and Discussion 


Structure Determination. To prepare NL63-CoV RBD for crys- 
tallization, we designed 36 RBD constructs with different N- and 
C-termini, on the basis of biochemical studies (19-21) and 
secondary structure prediction of the $1 sequence. We expressed 
and purified each of these fragments in insect cells. Alpha- 
chymotrypsin treatment of one of these fragments (residues 
461-616) yielded a smaller but more stable fragment (residues 
481-616), which was subsequently crystallized in complex with 
human ACE2 peptidase domain. We determined the structure 
by molecular replacement, using ACE2 as the search model 
(Figs. 1B and 2C). A 4-fold noncrystallographic symmetry 
averaging within the crystal and multiple cross-crystal averaging 
of the ACE2 region with 3 other structures (SARS-CoV-RBD- 
ACE2 complex, ACE2, and ACE2-inhibitor complex) (5, 24) 
dramatically improved the electron density in the NL63-CoV 
RBD region (Fig. 1C). We refined the structure at 3.3-A 
resolution (Table $1). The final model contains residues 19-614 
of human ACE2 and residues 482-602 (except for a disordered 
loop 555-565) of NL63-CoV RBD. The model also contains 
glycans N-linked to viral residues 486 and 512 and to ACE2 
residues 90 and 546. 


Structure of NL63-CoV RBD. The core of NL63-CoV RBD can be 
best described as a B-sandwich. It consists of 2 layers of 
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3-stranded B-sheets, one with mixed polarity (strands 3-1-2) and 
the other antiparallel (strands 4-5-6) (Figs. 1 B and D and 2 C 
and EF). The 2 layers stack tightly against each other through 
extensive hydrophobic interactions (Fig. 34). A short 3-stranded, 
antiparallel B-sheet (strands 2b-4b-6) stabilizes the distal end of 
the RBD (opposite the receptor-binding interface); this end 
contains both the N- and C-termini of the RBD and hence likely 
interacts with the rest of Sl and is not as ordered as the 
receptor-binding interface (Fig. 1B). Three disulfide bonds are 
present in the RBD (Fig. 34). The Cys-497—Cys-500 disulfide 
bond strengthens a critical receptor-binding loop (see below). 
Two other disulfide bonds, connecting Cys-516—Cys-567 and 
Cys-550—Cys-577, strengthen the region at the distal end. Ala- 
nine substitutions for Cys-516, Cys-550, Cys-567, or Cys-577 
decrease protein stability and abolish protein expression (Table 
S2) (21). Overall, the B-sandwich structure of NL63-CoV RBD 
is a unique protein fold; a search of the protein database using 
DALI (25) did not reveal any related structures. 

Three discontinuous receptor-binding sites on NL63-CoV 
RBD, which we term receptor-binding motifs (RBMs), are pre- 
sented by the B-sandwich core to bind ACE2. On the basis of 
their order in the primary structure, we define the 3 RBMs as 
RBM1 (residues 493-513), RBM2 (residues 531-541), and 
RBM3 (residues 585-590) (Figs. 1D and 3B). All 3 RBMs are 
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Fig.2. Structural comparison of NL63-CoV and SARS-CoV RBDs. (A) Domain 
structure of NL63-CoV S1. The boundaries of the RBD were determined by 
limited proteolysis of longer S1 fragments, followed by N-terminal sequencing 
and mass spectrometric analysis of the digestion fragments. The RBMs were 
identified from the crystal structure of the NL63-CoV RBD in complex with 
ACE2. (B) Domain structure of SARS-CoV $1 (5). (C) Another view of the 
structure of the NL63-CoV-RBD-ACE2 complex, which is derived by rotating 
the one in Fig. 1B by 90° clockwise along a vertical axis. (D) Structure of the 
SARS-CoV-RBD-ACE2 complex (PDB 2AJF) (5), from the same orientation as in 
C. (E) Schematic illustration of the topology of NL63-CoV RBD. Strands are 
drawn as arrows. (F) Schematic illustration of the topology of SARS-CoV RBD. 
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connected to B-strands and are thus B-loops. RBM3 is more 
compact than RBM1 and RBM2. Nevertheless, RBM1 is 
strengthened by the Cys-497—Cys-500 disulfide bond, and RBM2 
is reinforced by key hydrogen bonds involving Asp-538 (Fig. 3B). 
Alanine substitutions for Cys-497, Cys-500, or Asp-538 abolish 
receptor binding without significantly affecting protein expres- 
sion (Table $2) (21). At the top of the RBD, the 3 protruding 
RBMs surround a shallow bowl-shaped cavity, which is impor- 
tant for receptor binding (Fig. 44). 

NL63-CoV and SARS-CoV RBDs have no structural homol- 
ogy in either cores or RBMs, although they are both located in 
the S1 central regions (Fig. 2.4 and B). The core of SARS-CoV 
RBD is a single-layer, 5-stranded antiparallel B-sheet, stabilized 
by 3 short a-helices (Fig. 2 D and F); the core of NL63-CoV RBD 
is a B-sandwich with novel topology (Fig. 2 C and £). The 
number of B-sheet layers, the number of B-strands, and most 
notably, the topologies of the B-sheets all differ in the 2 RBD 
cores. The RBMs of SARS-CoV and NL63-CoV are also 
completely different. The SARS-CoV RBM is a continuous 
70-residue-long subdomain that lies on one edge of the core (Fig. 
2 D and F). The NL63-CoV RBMs are 3 discontinuous B-loops 
(Fig. 2 C and £). The lack of structural homology in RBD cores 
and, more importantly, the lack of structural homology in RBMs, 
indicate 2 independent ways in which NL63-CoV and SARS- 
CoV recognize their common receptor protein. 
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Fig. 3. Structural features of NL63-CoV RBD and RBMs. (A) Beta-sandwich 


core structure of NL63-CoV RBD. Hydrophobic residues between B-sheet layers 
are in yellow, cysteines in blue, glycans in green, and glycosylated asparagines 
in magenta. There exist 3 predicted N-linked glycosylation sites in the RBD 
(Asn-486, Asn-506, and Asn-512), and 2 of them are confirmed in the structure 
(Asn-486 and Asn-512). (B) NL63-CoV RBMs. (C) Residues on NL63-CoV RBMs 
that directly contact ACE2. 


Structure of NL63-CoV-Receptor Interface. NL63-CoV RBD binds 
to the outer surface of the N-terminal lobe of the ACE2 
peptidase domain (Figs. 1B and 2C). Three discontinuous virus- 
binding sites on ACE2, which we term virus-binding motifs 
(VBMs), are directly involved in virus binding (Fig. 4A and C). 
On the basis of their order in the primary structure, we define 
the 3 VBMs as VBM1 (residues 30-41 on the N-terminal helix), 
VBM2 (loop containing residues 321-330), and VBM3 (loop 
containing residues 353-356). The binding surfaces of ACE2 and 
RBD are complementary, with receptor VBMS3 inserting into the 


Fig. 4. Structural comparison of NL63-CoV-ACE2 and SARS-CoV-ACE2 in- 
terfaces. (A) Enlarged view of the NL63-CoV-ACE2 interface, from the same 
orientation as in Fig. 2C. VBMs on ACE2 are in blue, and RBMs on NL63-CoV are 
in red. Arrow indicates the bowl-shaped cavity surrounded by 3 viral RBMs. (B) 
Enlarged view of the SARS-CoV-ACE2 interface, from the same orientation as 
in A. (C) Footprint of NL63-CoV on the surface of ACE2. The view is derived 
from the one inA by rotating ACE2 by 90° along a horizontal axis, in such a way 
that the edge facing the viewer moves up. VBM1 residues are in orange, VBM2 
residues in magenta, and VBM3 residues in red. (D) Footprint of SARS-CoV on 
the surface of ACE2, from the same orientation as in C. VBM1b residues are in 
green. 
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Fig. 5. 


A common virus-binding hotspot on ACE2 for the binding of both 
NL63-CoV and SARS-CoV. (A) Stick-and-ball representation of the hotspot at 
the NL63-CoV-ACE2 interface. (B) Corey-Pauling-Koltun representation of the 
hotspot at the NL63-CoV-—ACE2 interface. (C) Stick-and-ball representation of 
the hotspot at the SARS-CoV-ACE2 interface. (D) Corey-Pauling-Koltun rep- 
resentation of the hotspot at the SARS-CoV-ACE2 interface. 


bowl-shaped cavity on RBD (Figs. 44 and 5A). All of the 
predicted and identified glycosylation sites on RBD and ACE2 
are located away from the virus-receptor interface (Fig. 34) (24) 
and thus do not interfere with virus—receptor interactions. A 
total of 11 viral residues directly contact 16 receptor residues 
(Figs. 3C and 4C and Table S3), and the binding buries approx- 
imately 1,300 A? total surface area. This binding interface is 
slightly smaller than other virus—receptor binding interfaces that 
typically range between 1,400 and 1,800 A? (2-6). Nevertheless, 
the binding affinity between NL63-CoV RBD and ACE2 is high. 
By surface plasmon resonance using Biacore, we measured the 
dissociation constant between the 2 proteins as 34.9 nM (Figs. 
1F and Fig. S1). 

NL63-CoV and SARS-CoV bind to common regions on 
ACE2, even though they have no structural homology in RBD 
cores or RBMs. Both viruses recognize the same 3 VBMs on the 
outer surface of the N-terminal lobe of ACE2 (SARS-CoV 
recognizes an additional short VBM1b) (Fig. 4 A and B). 
Compared with SARS-CoV RBD, NL63-CoV RBD has less 
extensive contact with VBM1 but more extensive contact with 
VBM2 and VBM3 (Fig. 4 C and D). In addition, nearly the same 
number of ACE2 residues directly contacts each of the 2 viruses. 
Among these ACE2 residues, 7 directly contact both of the 
viruses (Table S3). Overall, although SARS-CoV forms a larger 
binding interface with ACE2 than NL63-CoV does, the 2 viral 
RBDs bind to ACE2 with similar affinity (Figs. 1F and Fig. S1) 
(26). 


A Common Virus-Binding Hotspot on ACE2. A virus-binding hotspot 
on ACE2 lies at the center of the NL63-CoV-receptor interface. 
In the structure of unbound ACE2, Lys-353 projects into solution 
(24). Upon NL63-CoV binding, Lys-353 becomes embedded in 
a hydrophobic tunnel surrounded by 2 aromatic rings of ACE2 
Tyr-41 and NL63-CoV Tyr-498 and by 2 alkyl chains of ACE2 
Asp-37 and NL63-CoV Ser-535 (Fig. 5 A and B). At the end of 
the tunnel, ACE2 Asp-38 forms a salt bridge with Lys-353, 
neutralizing its charge. Because of the hydrophobic environ- 
ment, this salt bridge is energetically stabilizing (27) and critical 
for virus—receptor interactions. Alanine substitutions for Lys- 
353 or any other residues involved in the hotspot structure 
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abolish NL63-CoV binding (20, 21). The same virus-binding 
hotspot on ACE2 1s also key to the binding of SARS-CoV. The 
structures of the hotspots at the 2 different virus—receptor 
interfaces are strikingly similar (Fig. 5). The receptor parts are 
nearly identical, with subtle changes in protein side chain 
conformations. In the viral parts, Thr-487 and Tyr-491 on 
SARS-CoV replace Ser-535 and Tyr-498 on NL63-CoV, respec- 
tively, as 2 of the 4 tunnel walls. Adaptation to the hotspot on 
ACE2 is critical for SARS-CoV pathogenesis. A single T487S 
mutation that disturbs the hotspot decreases the binding affinity 
between SARS-CoV RBD and human ACE2 by more than 
30-fold and may explain why the SARS infections in 2003-2004 
were mild (23, 26, 28). Therefore, virus-binding hotspots on 
receptor proteins likely dictate virus—receptor interactions, viral 
pathogenesis, and viral transmissibility and thus are potentially 
major binding targets for viruses. 


Receptor-Recognition Mechanisms of Other Group | Coronaviruses. 
Despite their homologous spike proteins, NL63-CoV and other 
group I coronaviruses recognize different receptors. The S1 
central regions of porcine transmissible gastroenteritis virus 
(TGEV), feline infectious peritonitis virus (FIPV), porcine 
respiratory coronavirus (PRCV), and canine coronavirus 
(CCoV) are highly homologous to each other, and they all use 
APN from their respective host as receptors (29). On the basis 
of structural data for NL63-CoV RBD and sequence alignment 
of the S1 central domains of group I coronaviruses (Fig. S2), we 
make the following predictions. First, the B-sandwich core 
structure of NL63-CoV RBD 1s likely conserved among group I 
coronaviruses. The sequences of all of the B-strands are con- 
served, so are all of the 4 cysteines in the core that form 2 
disulfide bonds and strengthen the B-sandwich structure. Sec- 
ond, we predict that the 3 ACE2-binding RBMs in NL63-CoV 
are APN-binding sites in TGEV, FIPV, PRCV, and CCoV; 
structural changes in these RBMs result in the recognition of 
different receptors by these homologous viruses. Compared with 
NL63-CoV RBMs, the predicted RBM1 of other group I coro- 
naviruses has a 2-residue insertion and loses a critical disulfide 
bond, the predicted RBM2 has a 9-residue insertion and 2 
additional cysteines that likely form a disulfide bond, and the 
predicted RBM3 contains a 1-residue deletion. These structural 
changes are expected to abolish ACE2-binding capabilities of 
NL63-CoV and may help confer APN-binding capabilities of 
other group I coronaviruses. Structural studies of these RBDs 
and RBMs will further elucidate the receptor-recognition mech- 
anisms of these important animal coronaviruses. 


Evolution of Coronaviruses. Coronaviruses are believed to have 
common ancestors because they share similar replication mech- 
anisms, genomic structures, and overall gene sequences (10). 
Among all of the coronavirus genes, the one encoding the spike 
protein is the most variable. Between the spike protein subunits, 
S1 is more variable than $2. The current structural divergences 
of the S1 subunits reveal the tremendous evolutionary pressure 
that coronaviruses face to adapt to different host receptors, and 
they also reflect on the evolutionary history of coronaviruses and 
their receptor selections. 

Our study clarifies the baffling evolutionary relationships 
between NL63-CoV and SARS-CoV. First, because of the lack 
of structural homology in RBD cores and more importantly, the 
lack of structural homology in RBMs, NL63-CoV and SARS- 
CoV must have evolved their ACE2-binding capabilities inde- 
pendently, disproving the previous hypothesis that the newly 
emerged SARS-CoV acquired its RBM from the more ancient 
NL63-CoV or an NL63-CoV-related group I coronavirus (20, 
23). Second, binding to the same virus-binding hotspot on their 
common receptor protein, the 2 viruses present a remarkable 
case of functionally convergent evolution that occurred after 
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divergent evolution. Whereas differences in ancient host envi- 
ronments likely caused the 2 viruses to initially diverge, the 
virus-binding hotspot on their common receptor protein was the 
probable driving force for the subsequent convergent evolution. 


Materials and Methods 


Protein Purification and Crystallization. The NL63-CoV spike protein receptor- 
binding region (residues 48 1—616) containing an N-terminal honey bee melit- 
tin signal peptide and a C-terminal His tag was expressed in Sf9 insect cells 
using the Bac-to-Bac system (Invitrogen) and then purified as previously 
described (5, 30). In brief, the protein was harvested from Sf9 cell supernatants 
and loaded onto a Ni-NTA column. The protein was eluted from Ni-NTA 
column with imidazole and further purified by gel filtration chromatography 
on Superdex 200 (GE Healthcare). The protein was concentrated to 10 mg/mL. 
The peptidase domain of human ACE2 was expressed and purified using the 
same protocol as above. To purify the RBD-ACE2 complex, ACE2 was incu- 
bated with excess NL63-CoV RBD for 15 min at room temperature, and the 
complex was purified by gel filtration chromatography. 

Crystals of the RBD-ACE2 complex were grown in sitting drops at 12 °C, 
over wells containing 20% PEG6000 and 100 mM Na citrate pH 5.5. Crystals 
were harvested in 1 week, stabilized in 25% PEG6000, 100 mM Na citrate 
pH5.5, and 30% ethylene glycol, and flash frozen in liquid nitrogen. 


Structure Determination and Refinement. X-ray diffraction data were col- 
lected at Advanced Photon Source (Illinois) beamline 19 ID. The crystal 
contains 4 complexes per asymmetric unit. The structure was determined by 
molecular replacement using ACE2 peptidase domain structure as the 
search model (5). To overcome the model bias problem common to the 
molecular replacement method, we performed 4-fold noncrystallographic 
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symmetry averaging within the crystal and multiple cross-crystal averaging 
of the ACE2 region with 3 other structures (SARS-CoV-RBD-ACE2 complex, 
ACE2 by itself, and in complex with an inhibitor) (5, 24). This procedure 
dramatically improved electron density in the NL63-CoV RBD region. The 
structure of the NL63-CoV-RBD-ACE2 complex was refined to a final R¢ree of 
30.8% and Rwork of 27.6%. Data and refinement statistics are shown in 
Table S1. Software used for data processing, structure determination, and 
refinement is also listed in Table S1. 


Kinetics and Binding Affinity of NL63-CoV RBD and Human ACE2 by Surface 
Plasmon Resonance Using Biacore. The binding reactions between NL63-CoV 
RBD and human ACE2 were assayed by surface plasmon resonance using a 
Biacore 3000. ACE2 was immobilized on a C5 sensor chip directly. The surface 
of the sensor chip was first activated with N-hydroxysuccinimide; ACE2 was 
then injected and immobilized to the surface of the chip; last, the remain- 
ing activated surface of the chip was blocked with ethanolamine. Soluble 
NL63-CoV RBD was introduced at a flow rate of 20 uL/min at different 
concentrations. Kinetic parameters were determined with BlAevaluation soft- 
ware (Biacore) and are shown in Fig. 1E£. 
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